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1.  INTRODUCTION 


1.1  Purpose  of  the  Report 

This  report  describes  the  development  and  applications  of  a  computer  model  that  simulates 
steady  flow  inside  a  turboshaft  engine  test  facility.  This  model  was  developed  using  the 
Structured  Finite  Volume  (SFV)  code  [1,2],  a  general  purpose  computational  fluid 
dynamics  (CFD)  code.  The  report  provides  details  of  the  mathematical  model  and  reviews 
the  results  obtained  from  the  computations.  This  work  was  undertaken  by  Adaptive 
Research  Corporation,  for  the  Naval  Facilities  Engineering  Command,  under  Contract 
Number  N47408-91-C-1228. 


1.2  Background 

U.  S.  Navy  aircraft  engines  completing  maintenance  or  repair  are  operated  in  test  cells 
where  they  must  meet  performance  specifications  before  reinstallation.  Test  cells  are  fuUy 
enclosed,  sound  absorbing  hangers.  The  principle  function  of  these  facilities  is  to  provide 
repeatable,  specific  test  conditions.  In  the  process,  the  test  cells  must  attenuate  the  jet  and 
propeller  noise  sufficiently  for  personal  to  carry  out  normal  activities  near  the  facility,  cool 
and  control  the  engine  exhaust  flow  to  the  point  that  cell  components  are  not  damaged, 
and  meet  environmental  standards. 


The  Navy  has  recently  completed  construction  of  the  prototypes  of  a  new  generation  of  test 
ceils  designed  for  turboshaft  and  turboprop  engines.  These  test  cells,  located  at  the  Marine 
Corps  Air  Facility,  Camp  Pendleton,  California,  are  scheduled  to  begin  operating  shortly. 

1.3  Objective  of  the  Study 


The  objective  of  this  work  is  to  provide  a  computational  fluid  dynamics  model,  employing 
the  SFV  code,  of  these  prototype  engine  test  cells  for  use  in  performance  evaluation  and  for 
troubleshooting. 


iooesslon  For 

HTIS  GPASrI 


DTir  T.<iH 
Ur.r-!: -laMiioed 
j  J.i'jt ;  o.'il  1  oil-. 


DTIC  QUALITY  !N!]rE>'Tpr>  By _ _ 

Dl3%riVutl*n/ 

1  Availability  Codeo 

lAvall  and/or 
list  Speolel 


□  □ 


This  report  consists  of  six  main  sections.  Following  this  introduction,  Section  2  presents 
the  geometric  details  of  the  turboshaft  test  cell  along  with  operation  conditions.  Next, 
Section  3  contains  information  about  the  mathematical  basis  of  the  numerical  model.  In 
Section  4  the  numerical  details  of  the  model  such  as  boundary  conditions  are  presented. 
The  results  of  the  present  study  are  presented  in  Section  5  and  the  conclusions  drawn  are 
given  in  the  final  section. 


2.  DESCRIPTION  OF  THE  SIMULATED  PROCESS 


2.1  Geometry 

Details  of  the  geometry  are  provided  in  Appendix  A.  The  test  bay  is  rectangular.  The 
engine,  without  the  rotor,  is  mounted  on  a  cart  and  aligned  with  the  dynamometer  exhaust 
stack  and  the  augmenter  tube.  Also  on  the  cart  is  an  air  dynamometer  which  is  attached 
by  a  shaft  to  the  engine.  Air  for  the  engine  and  for  the  dynamometer  enters  the  cell 
through  two  rows  of  acoustic  vanes  located  in  the  front  wall  of  the  test  bay.  These  vanes 
are  not  arranged  symmetrically,  and  the  location  of  the  engine  does  not  fail  along  the 
center  line  of  the  test  cell.  Engine  exhaust,  along  with  some  cooler  test  bay  air  entrained 
by  the  jet,  leaves  the  facility  through  the  augmenter  tube,  discharging  to  a  chimney  and 
then  to  the  atmosphere.  The  dynamometer  has  its  own  exhaust  system.  Both  the  chimney 
and  the  dynamometer  exhaust  stack  are  provided  with  two  rows  of  acoustic  baffles. 

2.2  Operation  Conditions 

The  engine  considered  in  this  modeling  effort  is  the  T700  turboshaft  engine  manufactured 
by  General  Electric.  The  operation  conditions  for  the  engine  are  listed  in  Table  1. 

For  the  dynamometer,  a  mass  flow  of  14.0  Ib/s  (6.35  kg/s)  was  used.  The  average  exit 
temperature  used  for  the  dynamometer  was  413  F  (212  C).  The  ambient  conditions  outside 
the  cell  were  29.91  mm  Hg  (1.013  x  10^  N/m^l,  77  F  (25  C),  and  a  dry  bicomponenl 
atmosphere  (76.83  wt  per  cent  N2). 
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TABLE  1.  T700  Turboshaft  Engine  Operating  Conditions 


Exhaust  Flow 

11.2  Ib/s 

5.08  kg/s 

Fuel  Flow 

0.222  Ib/s 

0.101  kg/s 

Average  Exit  Temperature 

1200  F 

650  C 

Exhaust  Composition  (wt%) 

0.7512  N2 

0.1548  O2 

0.0690  CO2 

0.0250  H2O 
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3.  MATHEMATICAL  FORMULATION 


3.1  The  Structured  Finite  Volume  (SFy”)  Mathematical  Formulation 

The  SFV  computer  code  is  a  general-purpose  CFD  code  in  wide  use  throughout 
government,  industry,  and  academia.  It  can  model  two—  or  three-dimensional,  laminar  or 
turbulent,  single—  or  two— phase  fluid  flows  in  arbitrary  geometric  flow  domains.  SFV 
solves  the  conservation  equations  for  mass,  momentum,  and  energy.  The  details  of  this 
solution  procedure  follow. 

3.1.1  The  Governing  Equations 

The  flow  field  in  a  given  geometry  can  be  described  by  the  conservation  equations  for  mass, 
momentum,  and  energy.  These  equations  can  be  expressed  in  the  following  form 

1'!'+  V  (p  V  )  =  0  (i) 


+  V  (^W)  =  -  rp  +  T  W) 
^  +  7  (pVh)  =  75  +  +  WP 


where 


V 

is 

the 

p 

is 

the 

p 

is 

the 

^eff  ~ 

is 

the 

h 

is 

the 

q 

is 

the 

<P 

is 

the 

time-mean  velocity; 
gas  density; 
static  pressure; 
"effective"  viscosity;’ 

enthalpy; 

diffusive  energy  flux;  and 
dissipation  function. 


(2) 

(3) 
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The  effective  viscosity  is  defined  by  the  relation 


k 


where  /x  is  the  molecular  viscosity  and  the  turbulent  viscosity,  is  deduced  by  employing 

a  turbulence  model.  SFV  employs  the  two— equation  k— e  turbulence  model  [3].  This 
treatment  of  turbulence  requires  the  solution  of  two  additional  partial  differential  equations 


+  7(/)VK)  =  7 


^  =  7 


^eff 


76 


(4) 

(5) 


where 


k  -  is  the  turbulence  kinetic  energy; 

6  -  is  the  dissipation  rate  of  turbulence  kinetic  energy; 

Cp  C2,  model  constants;  and 

—  is  the  volumetric  production  rate  of  k,  defined  by 


Pi,  = 


dY.  dY.' 


dY. 

1 

W. 

J 


When  single— phase  three-dimensional  flow  is  under  consideration,  the  dependent  variables 
that  require  solution  are 


P 

u 

V 

w 

k,6 

h 


—  the  static  pressure; 

—  the  horizontal— direction  velocity  component; 

—  the  vertical— direction  velocity  component; 

—  the  axial-direction  velocity  component; 

—  turbulence  quantities  (described  above); 

—  the  static  enthalpy; 
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Cl 

C2 


—  the  mass  fraction  of  engine  gases;  and 

—  the  mass  fraction  of  dynamometer  gases. 


General  Form  of  the  Governing  Equations 


The  governing  equation  for  each  dependent  variable  can  be  reduced  to  a  single  general  form 
which  can  be  succinctly  represented  in  vector  notation  as 


difnp) 

“oTT 


(6) 


The  source  term  includes  both  sources  and  sinks  of  (fi  plus  any  other  terms  beyond  those 
which  appear  on  the  left-hand  side  of  the  equation. 


Because  of  the  generality  of  the  SFV  code,  a  standard  gradient-diffusion  law  is  provided 
for  the  diffusion  flux,  namely 


Where  diffusion  of  the  conserved  quantities  does  not  follow  the  standard  form  of  Equation 
7,  compensating  terms  are  added  to  the  source  term  S^.  In  that  equation,  is  the 

exchange  coefficient  for  the  variable  defined  as 


r 


'P 


^eff 

^^^,eff 


(8) 


where  the  effective  Prandtl/Schmidt  number  governing  diffusive  transport  of  the 

variable  ip. 


The  exception  to  the  general  form  given  by  Equation  6  is  the  continuity  (mass 
conservation)  equation.  The  pressure,  p,  has  been  classified  as  a  dependent  variable,  but 
does  not  appear  as  the  subject  of  a  transport  equation.  Instead,  the  pressure  is  associated 
with  the  continuity  equation  which  can  be  manipulated  to  derive  a  pressure— correction 
equation.  SFV  employs  a  variant  of  the  SIMPLE  pressure— correction  procedure  [4].  The 
transport  coefficient  and  source  term  are  provided  for  each  dependent  variable  in 

Table  2. 


TABLE  2.  Transport  Equations  and  Turbulence  Model  Constraints 


S 


Continuity  Momentum: 
X— direction 

y— direction 

z— direction 

Kinetic  Energy 
Dissipation  Rate 
Enthalpy 
Concentration 


1  0  0 


u 

^eff 

-  ^ 

v 

^eff 

-1 

w 

^eff 

-H 

k 

^t 

e 

h 

a 

e 

j(Cj  -Cj» 

h 

0 

C 

0 

C/i  =  0.09,  Cj  =  1.44,  Cj  =  1.92, 


=  1.0,  =  1.314 


3.1.3  General  Form  of  the  Finite-Domain  Equations 

Integration  of  the  generalized  conservation  equation  (ff)  yields  a  Dnite— difference  analog  for 
each  dependent  variable  for  each  control  cell  into  which  the  flow  geometry  is  divided. 
The  integration  process  is  carried  out,  and  the  terms  in  the  finite-domain  equations  are 
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assembled  in  the  manner  described  in  many  standard  texts  on  numerical  analysis  [5]. 

For  each  simulation  performed,  the  SFV  program  will  assemble  and  then  solve  a  set  of 
simultaneous  coupled  equations  of  the  general  form  (for  point  P) 

(9) 

where  the  subscripts  E,  W,  N,  S,  H,  and  L  represent  the  neighbor  points  in  space.  The 
coefficients,  a^^  and  others,  represent  the  effects  of  convection  and  diffusion.  is  the 

linearized  source  term.  The  central  coefficient  ap  is  defined  as 

ap— ap+a^Y+aj^+ag+ajj+a^^  (10) 

The  set  of  equations  (9)  are  assembled  and  solved  for  each  dependent  variable  in  sequence. 
The  solution  of  the  flow  field  then  proceeds  iteratively  to  steady— state. 

3.1.4  Auxiliarv  Relations 

To  complete  the  mathematical  formulation  of  a  problem  the  specification  of  additional 
information  is  required:  namely,  fluid  properties  and  boundary  conditions. 

SFV  allows  the  arbitrary  specification  of  fluid  properties.  The  models  resulting  from  this 
study  will  employ  fluid  properties  typical  of  the  conditions  present  in  the  Camp  Pendleton 
test  cells.  SFV  also  allows  several  different  types  of  flow  and  thermal  boundary  conditions 
These  will  be  appropriately  chosen  so  as  not  to  sacrifice  the  fidelity  of  the  models  and  also 
to  allow  for  easy  modification  in  the  course  of  the  parametric  studies. 
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4.  NUMERICAL  DETAILS  AND  PHYSICAL  MODELS 


4.1  Introduction 

This  section  of  the  report  describes  the  various  aspects  of  the  modeling  method  necessary 
to  simulate  the  flow  and  heat  transfer  inside  the  test  facility.  First,  the  computational  grid 
employed  in  this  simulation  is  described.  Next,  the  physical  properties  are  defined.  This  is 
followed  by  a  discussion  of  the  boundary  conditions  used  in  the  computational  model. 

4.2  Grid  Description 

The  first  step  in  setting  up  a  numerical  model  is  the  specification  of  the  computational 
grid.  For  this  study  a  three  dimensional  Body— Fitted  Coordinate  (BFC)  grid  was 
employed.  The  grid  for  this  modeling  effort  is  quite  complex.  In  order  to  simplify  the  grid 
generation,  a  procedure  has  been  developed  to  allow  for  relatively  easy  modifications  of  the 
test  cell  configuration. 

In  this  method  the  user  specifies  various  geometric  quantities.  The  standard  input  files 
then  create  a  set  of  data  files  for  the  EasyMesh  [6)  program  which  is  a  two  dimensional 
BFC  package.  The  user  must  then  run  this  grid  generation  package  to  create  a  set  of  two 
dimensional  grids.  The  standard  input  files  are  then  re-executed  to  produce  the  final  grid 
by  stacking,  blending,  and  rotation  of  the  two  dimensional  grid  planes. 

The  final  computational  grid  for  this  study  is  shown  in  Figure  1.  This  grid  has  40  cells 
across  the  test  bay,  35  cells  in  the  vertical  direction,  and  77  computational  cells  in  the  axial 
direction.  The  grid  is  clustered  in  areas  of  geometric  changes,  areas  of  expected  high 
shearing  action,  and  areas  of  heat  transfer.  In  Figures  2  through  4  outlines  of  the  internal 
components  are  presented. 

4.3  Properties  of  the  Fluid 

The  properties  of  the  gases  mixture  inside  the  test  cell  were  calculated  based  on  the  local 
conditions.  The  local  heat  capacity  was  computed  using  the  gas  mass  fraction  and  the  cell 
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temperature  using  JANNAF  [7]  polynomial  expressions.  The  gas  density  was  calculated 
using  the  ideal  gas  law  (requiring  local  pressure,  molecular  weight,  and  temperature).  The 
laminar  kinematic  viscosity  was  assumed  to  have  a  constant  value  of  1.08  x  10  ft“/s  (1.0 
X  10~^  m^/s).  The  turbulent  viscosity,  which  tends  to  dominate,  was  calculated  from  the 
local  turbulence  quantities. 

Heat  transfer  through  the  augmenter  tube  wall  required  the  specification  of  the  material 
thermal  conductivity.  Inside  the  room  the  thermal  conductivity  of  mineral  fiber  with  a 
value  of  0.22  Btu/ft-s-R  (0.38  W/m-K),  was  used.  In  the  chimney  the  conductivity  of 
steel  was  used.  The  value  of  steel  was  taken  to  be  26  Btu/ft— s— R  (45  W/m— K). 

4.4  Boundary  Conditions 

For  this  model,  the  boundary  conditions  requiring  specifications  were 

•  the  engine, 

•  the  dynamometer, 

•  inlet  into  test  cell, 

•  outlet  of  dynamometer  stack  and  chimney, 

•  wall  boundaries,  and 

•  heat  transfer  through  the  augmenter  wall. 

These  boundary  conditions  and  sources  are  discussed  in  the  following  sections. 

4.4.1  Engine 

The  engine  is  modeled  as  a  hollow  cylindrical  shape  with  a  variable  cross  sectional  area. 
Inside  this  region  a  plate  of  computational  faces  were  blocked.  On  the  intake  side  a  mass 
sink  is  applied  which  corresponded  to  the  exit  mass  Qow  minus  the  fuel  intake. 

On  the  exhaust  side  a  mass  source  was  applied  to  account  for  the  engine  mass  flow.  Other 
boundary  conditions  applied  the  appropriate  sources  for  the  momentum,  enthalpy, 
turbulence,  and  concentrations.  The  exit  turbulence  quantities  were  calculated  assuming  a 


15  per  cent  turbulence  intensity  from  the  following  equation: 

k  =  1/2  (vij) 

e  =  0.1643  k*-^/l 
'  m 


where 

—  is  the  turbulence  intensity  and 

Ij^  —  is  the  mixing  length  (based  on  a  gap  distance  or  inlet  radius) 

4.4.2  Dynamometer 

Similarly  for  the  dynamometer,  a  plate  inside  the  dynamometer  was  blocked.  On  the 
inside  face  a  mass  sink  was  applied.  A  mass  source  was  applied  to  the  top  face  along  with 
sources  for  the  momentum,  enthalpy,  turbulence  and  concentrations.  A  turbulence 
intensity  of  10  per  cent  was  assumed. 

4.4.3  Inlet 

.A.t  the  inlet  a  fixed  pressure  was  applied  which  corresponded  to  the  ambient  pressure. 

From  this  and  the  other  boundary  conditions  the  mass  flow  was  computed.  Since  this  mass 
flow  varies  during  the  solution  process,  the  momentum  source  was  calculated  from  the 
adjusting  inflow.  The  enthalpy  and  concentrations  were  considered  to  be  ambient.  A 
turbulence  intensity  of  2  per  cent  was  used  to  calculated  the  incoming  turbulence  values. 

A  constant  velocity  of  3.3  ft/s  (1.0  m/s)  was  used  in  the  turbulence  calculations. 

The  bafQes  located  at  the  inlet,  dynameter  exhaust,  and  chimney  were  not  directly 
modeled.  These  effects  were  modeled  by  calculating  the  momentum  loss  through  the 
baffles  from 


where 


Ap  =  0.5  K/jv^ 

K  is  a  loss  coefficient. 
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Supplied  with  pressure  drop  across  the  baffles  along  with  gas  density  and  velocity  an 
accurate  determination  of  K  is  possible.  Pressure  data  only  existed  for  the  baffles  of  the 
exhaust  stack.  However,  this  value  was  for  a  location  inside  the  baffles.  It  was  assumed 
that  the  pressure  drop  across  the  baffles  was  double  this  value.  Using  estimates  for  density 
and  velocity  a  K— loss  factor  of  0.2  was  calculated.  This  was  lowered  to  0.1  for  ail  baffles  in 
an  effect  to  better  predict  temperature  data  (Section  5.4). 

4.4.4  Outlet 

At  the  two  ’outlets’  a  fixed  pressure  was  used.  If  no  recirculation  occurs  this  is  all  that  is 
required  as  in  the  case  of  the  dynamometer  stack.  However,  in  the  chimney,  the  possibility 
of  recirculation  at  the  exit  was  extremely  great.  Momentum  and  energy  sources  were 
supplied  if  and  when  recirculation  at  the  exit  plane  was  calculated.  The  momentum  source 
was  based  on  the  amount  of  mass  being  brought  into  the  domain  and  the  energy  source  was- 
based  on  ambient  gas  conditions. 

4.4.5  Wall  Boundaries 

At  the  wail  boundaries  the  wall  function  approach  [3]  was  used  to  account  for  momentum 
losses.  This  condition  is  not  applied  directly  to  the  wall  but  rather  at  a  point  outside  the 
viscous  sublayer,  where  the  logarithmic  law  of  wall  prevails.  This  is  where  turbulence  is 
assumed  to  be  in  local  equilibrium  (i.e.,  the  generation  and  dissipation  of  turbulence  energy 
are  equal).  Heat  transfer  was  not  considered  except  through  the  augmenter  tube  wall. 

4.4.6  Heat  Transfer  through  Augmenter  Tube 

To  compute  the  heat  transfer  through  the  augmenter  tube,  a  composite  heat  transfer 
treatment  was  used.  Surface  heat  transfer  coefficients  on  both  surfaces  (inside  and  outside) 
were  calculated  using  the  standard  SFV  wall  fiiction  treatment.  The  overall  heat  transfer 
conductance  between  the  fluid  cell  adjacent  to  the  inner  wall  and  the  one  adjacent  to  the 
outer  wall  was  evaluated  using  these  two  surface  coefficients  (in  the  form  of  a  heat  transfer 
resistance)  and  the  assumed  thermal  resistance  through  the  tube  wall.  This  last  resistance 
was  evaluated  using  a  conductivity  of  insulation  and  thickness  of  insulation  provided  by 
the  Navy.  In  the  chimney,  where  there  was  no  insulation,  the  conductivity  and  thickness 
of  the  augmenter  tube  wall  was  used  for  the  solid  contribution  to  the  overall  heat  transfer 
conductions. 
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5.  DISCUSSION  OF  RESULTS 


This  report  deals  mainly  with  results  from  one  computational  run.  Differences  in  the 
temperatures  predicted  and  measured  values  gave  rise  to  the  investigation  of  various  items 
in  regard  to  their  effect  on  temperature.  These  items  will  also  be  discussed  in  the  following 
sections. 

5.1  Velocity 

The  velocity  vectors  are  presented  in  Figures  5  through  10.  A  cross  sectional  view  of  the 
vectors  are  show  in  Figure  5.  This  plot  provides  a  good  overall  display  of  the  major  flow 
patterns  within  the  cadculation  domain.  Two  plots  of  velocity  around  the  dynamometer  are 
shown  in  Figures  6  and  7.  The  first  one  again  shows  a  side  view.  In  this  plat  the  major 
features  of  the  dynamometer  and  engine  inlet  are  clearly  visible.  This  plot  shows  the  air 
movement  into  the  two  10"  diameter  (0.254  m)  openings  of  the  dynamometer  and  also  the 
14"  diameter  (  0.366  m)  engine  inlet.  This  plot  shows  the  alignment  of  the  dynamometer 
with  the  dynamometer  exhaust.  In  this  case  the  dynameter  is  shifted  1"  (0.0254  m) 
forward.  With  this  alignment  part  of  the  flowing  exiting  the  dynamometer  actually  enters 
into  the  dynamometer  exhaust  sleeve.  In  Figure  7  another  section  (X— Y)  across  the 
dynamometer  is  presented.  In  this  plot  note  the  line  above  the  circular  inlet.  This  is  the 
location  of  the  boundary  conditions  to  account  for  the  increase  in  temperature  in  the 
dynamometer.  Because  of  the  swirl  in  the  dynamometer,  the  exit  velocity  components 
were  supplied  at  a  45  degree  angle.  The  velocity  for  the  X  and  Y— direction  was  59.1  ft/s 
(18.0  m/s). 

Figure  8  shows  the  velocity  vectors  around  the  gap  between  the  engine  and  augmenter 
tube.  The  entrainment  of  the  surrounding  gases  into  the  tube  is  clearly  shown.  The 
maximum  velocity  of  the  engine  is  438  ft /s  (133  m/s).  The  velocity  in  the  chimney 
section  is  shown  in  Figure  9.  In  the  wall  region  the  ducting  changes  from  a  circular  duct 
into  a  square  duct.  The  axial  momentum  propels  the  gases  inside  the  augmenter  tube  to 
the  aft  side  of  the  bend  section.  The  gases  do  not  have  enough  time  to  spread  to  the  entire 
cross  section  of  the  chimney,  thus  causing  a  recirculation  zone  at  the  exit  of  the  calculation. 

The  final  vector  plot  (Figure  10)  shows  a  top  side  view  of  the  velocity  vectors  entering  into 

the  test  cell.  The  entrance  velocity  is  approximately  5  ft/s  (1.5  m/s).  The  spreading  of 

the  incoming  air  and  the  recirculation  zone  beside  the  inlet  baffies  can  be  seen. 
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5.2  Temperature 

Temperature  plots  are  provided  in  Figure  11  through  19.  The  overall  temperature  field  is 
shown  in  the  Figure  11.  The  location  and  effect  of  the  two  heat  sources  can  be  observed. 
Details  of  the  temperature  field  around  the  dynamometer  are  displayed  in  Figures  12  and 
13.  The  exit  temperature  of  the  d3rnamometer  used  in  this  study  was  413  F  (212  C).  In 
Figure  12  the  effect  of  the  misalignment  between  the  dynamometer  and  the  dynamometer 
sleeve  can  be  noted.  If  not  for  the  skirt,  the  gases  would  enter  into  the  test  cell  with  a 
velocity  approximately  the  same  magnitude  as  leaving  the  dynamometer.  Originally  no 
swirl  was  used  as  a  boundary  condition.  Using  this  boundary  condition  no  temperature 
gradient  was  predicted  at  the  exit  (i.e.,  uniform  413  F  exit  temperature).  During 
additional  studies  described  in  Section  5.4,  a  decision  was  made  to  implement  a  swirling 
component  of  velocity. 

The  skirt  channels  the  flow  downward  back  into  the  suction  of  the  dynamometer.  The 
effect  of  the  swirl  at  the  dynamometer  boundary  condition  can  be  seen  in  Figure  13.  With 
this  boundary  condition  flow  is  entrained  on  the  right  side  (looking  fore  to  aft)  and  expelled 
on  the  left  side.  The  effect  of  the  entrainment  can  be  seen  in  the  temperature  plot  (Figure 
14)  at  the  exit  of  the  chimney  exhaust  section.  This  temperature  range  at  the  exit  varies 
from  386  F  to  413  F  (196  C  to  212  C). 

The  temperature  around  the  e.xit  of  the  engine  is  shown  in  Figure  15.  The  e.xit 
temperature  of  the  engine  was  1200  F  (650  C).  This  temperature  is  then  reduced  as  the 
engine  gases  are  mixed  with  the  entrained  air  in  the  augmenter  tube  as  shown  in  Figure  16. 
The  plot  of  temperature  at  the  exit  of  the  augmenter  tube  is  given  in  Figure  17.  At  this 
point  the  temperatures  range  from  681  F  to  802  F  (361  C  to  428  C).  The  average 
temperature  is  approximately  763  F  (406  C).  In  Figure  18  the  temperatures  at  the  exit  of 
the  chimney  section  are  plotted.  The  zone  of  recirculation  is  clearly  noted  by  a  region 
where  the  temperature  is  77F  (25  C),  which  corresponds  to  ambient  conditions.  An  axial 
temperature  profile  plot  is  provided  in  Figure  19.  In  this  plot  the  temperature  at  three 
locations  falling  inside  the  augmenter  tube  is  shown.  The  near  wall  temperatures  are 
shown  by  lines  (Y— cell  index)  labeled  12  and  24  while  the  core  flow  of  the  tube  is 
represented  by  line  19. 
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5.3  Other  Quantities 

In  Figure  20  the  pressure  field  located  approximately  2  ft  (0.6  m)  from  the  right  side  is 
displayed.  There  is  a  cell  depression  through  out  this  entire  computational  plane,  however 
in  terms  of  inches  of  water  the  values  are  quite  small.  Figure  21  shows  the  turbulence 
kinetic  energy  (TKE)  in  the  gap  region  of  the  engine.  Using  a  turbulence  intensity  value  of 
15  per  cent  an  inlet  value  for  TKE  of  202  m  /s  was  obtained. 

Two  global  numbers  for  the  calculation  were  calculated.  These  were  the  pumping  ratios  for 
the  dynamometer  and  the  engine.  The  values  calculated  for  this  computational  run  were 
-0.067  and  0.755,  respectively.  The  negative  value  for  the  dynamometer  pumping  ratio 
indicates  that  more  leakage  than  entrainment  occurred. 

Heat  transfer  through  the  augmenter  tube  was  at  the  rate  of  1.2  BTU/s  (1270  J/s). 
Approximately  90  per  cent  of  this  occurred  in  the  uninsulated  section  in  the  chimney. 

5.4  Comparison  With  Data 

During  October  of  1991  experimental  data  was  taken  for  the  test  case  that  was  simulated. 
The  results  are  listed  in  Table  3.  Included  in  the  table  are  the  values  obtained  from  the 
computation  run. 

T.ABLE  3.  Summary  of  Experimental  Results 


EXP 

CFD 

Pumping  ratio  engine 

0.68 

0.76 

Pumping  ration  dynamometer 

0.74 

-0.07 

Cell  depression  (in  H20) 

0.15 

<0.01 

Engine  augmenter  exit  average  temperature  (F) 

618 

763 

Engine  augmenter  exit  maximum  temperature  (F) 

628 

802 

Engine  exhaust  stack  average  temperature  (F) 

555 

nc* 

Engine  exhaust  stack  maximum  temperature  (F) 

596 

744 

Dynamometer  exhaust  duct  average  temperature  (F) 

348 

<400 

Dynamometer  exhaust  duct  maximum  temperature  (F) 

362 

413 

Dynamometer  exhaust  exit  maximum  temperature  (F) 

*Not  calculated 

413 

413** 

**  Boundary  condition 
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A  greater  degree  of  confidence  was  given  by  NCEL  to  the  actual  temperature 
measurements  than  other  quantities  (i.e.,  pumping  ratios,  velocity  measurements).  The 
model  seemed  to  overpredict  the  values  of  the  temperatures  for  both  the  dynamometer  and 
augmenter  exists.  Because  of  this  several  weeks  were  spent  investigating  the  differences. 

The  dynamometer  was  the  first  area  to  be  studied.  Given  the  geometry  with  a  gap  width 
of  only  1.25”  (0.0318  m)  and  the  misalignment  of  the  exhaust  and  the  stack,  it  was  difficult 
to  explain  this  amount  of  entrainment  indicated  by  the  data.  In  order  to  study  this  region 
a  smaller  3— dimensional  model  was  extracted  &om  the  larger  model.  During  this  study 
eleven  computational  tests  were  conducted  to  study  the  effect  of  the  turbulence  intensity, 
K~loss  factor,  exit  velocity  direction,  and  gap  space. 

Initially,  no  swirl  was  assumed  as  a  boundary  condition  for  the  dynamometer.  With  this 
condition,  no  entrainment  wai  observed  (only  a  leakage  into  the  test  cell).  Variations  of 
the  turbulence  intensity  and  the  K— loss  factor  produced  the  similar  results  (i.e.,,  no 
temperature  gradient  at  the  exit).  The  value  for  the  K^oss  factor  had  a  large  effect  on  the 
amount  of  leakage  between  the  gap.  The  larger  the  K-loss  factor  the  greater  the  leakage. 

It  was  only  when  swirl  was  added  to  the  dynamometer  boundary  condition  that 
entrainment  was  observed.  However,  in  all  cases  studied  there  was  a  negative  pumping 
ratio,  indicating  more  leakage  than  entrainment.  The  gap  between  the  dynamometer  exist 
and  the  exhaust  was  increased  over  several  cases  until  it  fell  0.5"  (0.0127  m)  below  the 
skirt.  The  major  result  was  that  as  the  gap  increase  so  did  the  amount  of  leakage  into  the 
test  cell. 

Under  no  condition  studied  did  the  model  match  the  data.  The  average  exit  temperature 
was  always  higher.  These  findings  point  to  explanations  or  combinations  as  to  the 
difference.  The  first  is,  that  actual  geometry  may  not  have  been  simulated  (i.e.,,  the  gap 
was  larger  or  a  different  misalignment).  The  second  is  the  boundary  condition  used  in  the 
dynamometer.  In  order  to  rigorously  model  the  dynamometer,  momentum  sources  would 
have  been  supplied.  With  the  proper  sources  the  model  would  have  accurately  predicted 
the  mass  flow  through  the  dynamometer.  However,  this  would  have  been  less  stable  and 
difficult  to  match  the  known  mass  flow  through  the  dynamometer.  Due  to  these  reasons  a 
plate  location  was  assumed  inside  the  dynamometer  (described  in  Section  4.4.2).  On  one 
side  of  the  plate  a  mass  sink  was  applied,  while  on  the  top  side  a  mass  source  was  applied 
which  took  into  account  the  temperature  rise  inside  the  dynamometer  If  this  boundary 
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condition  greatly  overpredicts  the  pressure  in  the  gap  region,  an  under  prediction  of  the 
entrained  mass  will  occur.  This  can  be  further  investigated  with  a  small  3— dimensional 
model. 

The  exit  temperatures  of  the  augmenter  tube  were  then  studied  in  detail.  During  these 
studies  the  following  four  parameters  were  examined;  engine  turbulence  intensity,  grid 
resolution,  K-4oss  of  inlet  and  chimney,  and  engine  temperature.  The  preliminary  results 
indicated  a  need  for  a  higher  engine  pumping  ration  to  account  for  an  over  prediction  in 
temperatures.  The  first  test  involved  lowering  of  the  engine  turbulence  intensity.  The 
result  to  the  pumping  ration  was  minimal. 

Grid  resolution  was  the  next  item  investigated.  A  smaller  3— dimensional  area  around  the 
augmenter  tube  was  assembled  from  the  larger  model.  Two  cases  were  conducted.  The 
first  used  the  same  grid  spacing  while  the  second  increased  the  number  of  cells  in  the  gap 
region  in  all  three  directions.  Preliminary  results  showed  little  difference  between  the 
pumping  ratios,  indicating  a  satisfactory  level  of  grid  refinement. 

At  this  point,  the  full  model  was  restarted  with  modifications  to  the  boundary  conditions 
namely;  swirl  at  the  dynamometer,  turbulence  intensity  of  engine  increased  from  8  to  15 
per  cent,  and  a  lower  K-4oss  factor  for  the  inlet  and  chimney  (note  these  are  the  results 
presented  in  this  report).  It  was  reasoned  that  by  lowering  the  K— loss  factor,  more 
entrainment  would  occur  as  in  the  case  of  the  dynamometer.  However,  it  was  found  that 
the  model  is  fairly  insensitive  to  the  K— ioss  factors  used  for  the  inlet  and  the  chimney  (i.e.. 
very  little  difference  in  the  pumping  ratio).  Reducing  the  K— loss  factor  from  earlier 
estimates  to  a  value  of  0.1  resulted  in  the  reduction  of  the  ceil  depression  from  0.008  inches 
of  H2O  to  0.004  inches  of  H2O. 

It  was  originally  thought  that  the  value  of  the  engine  turbulence  intensity  had  little  effect 
on  the  results,  as  indicated  by  the  first  test.  A  2— dimensional  model  was  developed  by 
NCEL  in  order  to  investigate  the  temperatures  within  the  augmenter  tubes.  This  model 
showed  that  the  mixing  inside  the  augmenter  tube  was  greatly  affected  by  engine 
turbulence  intensity,  as  shown  in  Figures  22  through  24.  In  these  cases  the  inlet  value  of 
turbulence  kinetic  energy  was  increased  from  65  to  520  to  977  m  /s  .  The  results  show 
that  as  the  intensity  is  increased  the  overall  mixing  inside  the  augmenter  tube  is  enhanced. 
However,  the  pumping  ratios  only  varied  between  0.428  and  0.463  for  the  three  cases. 
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closer  look  at  the  data  from  the  first  test  confirmed  the  effect  of  turbulence  intensity  on 
temperature  gradient. 

The  only  parameter  left  to  explain  the  differences  in  temperatures  was  the  value  used  for 
the  engine  exhaust  temperature.  A  valve  of  1200  F  (650  C)  had  been  used.  This  value  was 
obtained  from  engine  specifications  and  was  not  a  measured  quantity.  At  this  time  a 
1— dimensional  calculation  was  performed.  In  this  study  two  streams  were  "mixed" 
simulating  the  mixing  of  engines  gases  and  ambient  air  with  a  pumping  ratio  of  0.760.  In 
order  to  obtain  a  mixture  with  a  average  temperature  of  620  F  (327  C),  as  indicated  by  the 
data,  the  exit  engine  temperature  had  to  be  970  F  (521  C)  or  the  pumping  ratio  had  to  be 
increased  to  1.28  (at  1200  F).  The  large  model  was  repeated  with  the  engine  temperature 
lowered  to  1030  F  (554  C).  Figure  25  shows  the  plot  of  temperature  at  the  exit  of  the 
augmenter  tube. 

With  these  studies  of  the  augmenter  tube,  there  is  no  parameter  that  can  be  adjusted  to 
increase  the  pumping  ratio  from  0.760  to  1.28.  The  best  possible  explanation  for  the 
difference  in  temperature  is  that  the  exit  engine  temperature  is  actually  lower  than  the 
value  used.  If  the  profile  is  not  matched,  these  tests  indicate  that  the  value  of  the 
turbulence  intensity  can  be  modified  to  inaease  or  decrease  the  gradient  across  the  exit  of 
the  augmenter  tube.  If  the  locations  of  the  temperature  probes  were  to  be  laid  out  over  the 
predicted  values  from  the  model  (Figure  17),  it  can  be  noted  that  they  would  generally  fall 
in  an  area  of  average  temperature  of  the  plane.  The  probes  are  not  placed  in  a  location 
that  will  pick  up  the  full  gradient  of  the  exit  plane. 

The  measured  and  predicted  results  are  also  different  in  one  other  location.  Measured 
velocities  show  that  in  the  chimney  section  the  flow  exits  the  fore  side.  The  computational 
model  predicts  that  the  flow  exits  the  aft  side  with  recirculation  occurring  on  the  fore  side. 
Previous  turn  around  duct  studies  [8]  predict  (and  match  data)  the  transfer  of  momentum 
to  the  aft  side.  There  is  no  explanation  for  the  difference  in  results. 
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6.  CONCLUSION 

A  3— dimensional  model  of  a  Naval  test  cell  facility  has  been  developed.  This  model  allows 
the  user  to  investigate  various  parameters  of  interest  including  variations  in  geometry, 
ambient  conditions,  and  boundary  conditions.  The  results  of  this  study  show 

•  an  engine  pumping  ratio  of  0.755; 

•  an  dynamometer  pumping  ratio  of  -0.067; 

•  cell  depression  of  less  than  0.01”  of  water; 

•  an  augmenter  tube  average  temperature  of  763  F;  and 

•  an  exhaust  stack  temperature  near  the  value  of  the  dynamometer  exit. 

.Although  the  predicted  temperatures  were  higher  than  measured  values  possible 
explanations  to  the  differences  include 

•  difference  in  modeled  and  actual  geometries  and 

•  types  and  values  of  boundary  conditions  including  an  unmeasured  engine  exit 
temperature. 

Results  of  this  model  were  sensitive  to 

•  K— loss  value  for  the  dynamometer  exhaust, 

•  Swirl  at  augmenter  e.xit, 

•  Engine  turbulence  intensity,  and 

•  Engine  exit  temperature. 

Results  of  the  model  were  not  (or  least)  sensitive  to 

•  K-loss  for  inlet  and  chimney, 

•  Grid  refinement,  and 

•  Heat  loss  through  the  augmenter  tube. 
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Velocity  Vectors  Around  Gap  Region 
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Figure  10.  Velocity  Vectors  Inlet  Section 
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Figure  14.  Temperature  Exit  of  Dynamometer  Section 
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Temperature  Around  Gap  Region 
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